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Abstract 

By using similarity transformations approach, the exact propagator 
for a generalized one-dimensional Fokker-Planck equation, with linear 
drift force and space-time dependent diffusion coefficient, is obtained. 
The method is simple and enables us to recover and generalize special 
cases studied through the Lie algebraic approach and the Green function 
technique. 

The Fokker-Planck (FP) equation is a very useful tool in modelling a large 
number of stochastic phenomena in non equilibrium statistical physics as well 
as in chemistry, biology and many other disciplines of science [jl], ^, [|. The 
general one dimensional FP equation is given by 

J^P (x, t) = - (x, t) P (x, t) + /i 2 (X, t) P (X, t) . (1) 

In this work we are interested in studying the case where fi 1 and /i 2 (x,t) 
are linear functions of x with time dependent coefficients; namely 

Hi(x,t) =C(i)x + D(t), and /i 2 (x,t) = B(t) + E{t)x. (2) 

Here we assume that the coefficients B(t), C(t), E{t) are arbitrary real functions 
of time and D(t) satisfies some constraint that will be stated later so as to get 
exact analytic solutions to the problem. Furthermore, since /i 2 (x, t) must be 
positive semi-definite, the solution has to be limited either to the half axis 
x > for E(t) > or to the half axis x < for E(t) < 0. 

By these restrictions (2), the FP equation can be written 

^P(x,t)=H(x,t)P(x,t), (3) 
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where the operator H{x,t) (hereafter called Hamiltonian) is given by 



H(x, t) = -C{t) - [D{t) - 2E(t) + C(t)x] — + [B(t) + E(t)x] —. (4) 

Note that Eq. (||) includes two special cases: the FP equation with a linear drift 
force for E(t) = which was solved in the framework of Lie algebraic approach 
[|] and by the Green function technique ||, and the generalized FP equation 
with linear drift force for B(t) = and D(t) = ^E(t) {§. 

If the coefficients B,C,D, and E are constant, the solution of Eq.(||) may be 
given by the standard method using Laplace of Fourier transforms of P(x, t) with 
respect to x. This leads to a new function, say P(k, t), that obeys a first order 
partial differential equation and can be solved by the method of characteristics. 

If these coefficients are time dependent function, the solution of Eq. (||) may 
be written in the form of an ordered integral on time as 

P(x,t) = U(x,t)P(x,0), (5) 

where U{x 1 t) is the evolution operator. Since, H(x,t) does not commute with 
H(x,t') for t ^ t', U(x,t) reads 



U{x,t) = Texpj^ <ItH(x,t)^, 



with T being the time-ordering operator. To get over this difficulty, the familiar 
approaches use the Lie algebra j|, [| or the theorem of Suzuki [Q . 

In this paper, we propose to give the solution of this equation via the 
similarity transformations approach ^ ^| which consists in performing a series 
of transformations on Eq. (^) in order to write the evolution operator as a 
product of non time-ordering operators whose their action on P (x, 0) is known. 

Performing the transformation P(x,t) = X(x,t) Q(x,t) , where X(x,t) is 
an arbitrary invertible operator, one gets the new equivalent equation 

-Q(x,t)=H'(x,t)Q(x,t), (6) 
with the new Hamiltonian 

d 

H'(x,t) = X-\x 1 t)H(x,t)X(x 1 t) - X- 1 (x,t)—X(x,t). (7) 

Notice that since our so-called Hamiltonian is not necessarily Hermitian as in 
the quantum mechanical Schroedinger equation, the unitary condition on the 
operator X (X* = X^ 1 ) is not necessary. The operator X may have just to be 
invertible, in this case the transformations are called similarity transformations. 

Our first task is to eliminate the x-j^ operator in H(x,t). Thus, we choose 
X(x, t) as 

X(x,t) = exp(-a(t))exp f-a(t)x-^-\ with a(t) = J <2t(7(t), (8) 
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such that H(x, 0) = H'(x, 0) and P(x, 0) = Q(x, 0). 
Since 



i 9 v - „-a(*)A 



X" 1 a;X = e a w^, X~ 1 —X = c 

dx 



dx' 



it is easy to show that 



H'(x, t) = -e" Q W \D(t) - 2E(t)\ i- + e~ a(t) e~ a ^B(t) + E(t)x 

ox L 



dx 2 



(9) 



(10) 



When £7(4) = 0, and B(t) > 0, H'{x,t) reduces to a sum of commuting 
operators and -§^s)- The equation (g) is easily integrated to give 
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and consequently 



dx 







Q(x,t)=expi-d(t)—}e X pn(t)-^}Q(x,0), 



d 2 



dx 



P(x,t) = exp(-a(i)) exp -a(4)z— exp -d(4) — exp b(t)— P(x,0) 



dx 2 







(11) 



dx 



dx 2 J 



where 



6(4) 



/ dTB{T)e- 2a(T \ d(t) = [ drD(T)e 
Jo Jo 



■o(t) 



Using the well-known relations 

e^*/(x) = /(x + 7 ), e^/(z) = /(eV); V 7 . 
as well as the formula 



1 



/47T7 J _ 



dy exp 



1 

47 



y 



f(Vh) 



valid for Re 7 > and e = zbl, Eq.([l2|) may be put into the form 



P{x,t) 



dyK(x,fr,y,0)P(y,0), 



(12) 
(13) 

(14) 
(15) 

(16) 



where the transition probability density (propagator) K(x,t;y,0) is given by 



K(x,t;y,0) = 



exp 



(xe~ a ^ - d{t) - y y 



a (4) 



(17) 



This result agrees with the one obtained in refs ||, . Notice that since this is 
a Gaussian distribution, it may be easily checked by evaluating the average (x) 
and the moment (a; 2 ) directly from Eq. ([[]). 
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Let us now consider the case where E(t) ^ 0. In order to simplify Eq. (ph 

ft 2 t 
let us transform it by eliminating the operator in the Hamiltonian H by 

means of the unitary transformation defined as 

Q{x, t) = exp ( a{t)^- ] R(x, t), (18) 



dx 

with. a(t) — -§^ye~ Q (*). It is easy to show that the Hamiltonian H' becomes 

H" = -e~ a W (D(t) - 2E(t)+ a {t)e a{t) ) ^- + e~ a ^E(t)x-^. (19) 
V J ox ox z 

If we choose D(t) as 



D(t) = ^E(t)-a(t)e a V, (20) 



Eq. (W% is reduced to 



H » = e -«<.» E(t) (i|. +I ^). (21) 

The transformed equation 

^- t R(x,t)=H"(x,t)R(x,t), (22) 

is also a Fokker-Planck equation with a diffusion coefficient e~ a ^ E(t)x, which 
must be positive semi-definite for all x and t. Thus the solution of (|22| ) must 
be restricted to the half axis x > for E(t) > or to the half axis x < for 
E(t) < 0. This means that for any time there is no chance for the process to 
reach the negative half axis for E(t) > or the positive half axis for E(t) < 0. 
Thus R{x, t) is such that 

R(x, t) = for x < if E(t) < 0, or for x > if E(t) > 0. (23) 

Setting E(t) = e \E(t)\ with e = ±1 being the sign of E(t),the formal solution 
may be put as 



R(x, t) — 6{ex) exp 



.19 d 2 
e{t ^2dx- +X W 



R(x,0), (24) 



where e(t) = J* d T E{ T )e- a ^ . 

However, the usual method for solving equations of type ( p2] ) with the con- 
dition ( p3| ) is by the use of Green's function technique, via Laplace of Fourier 
transforms, and the method of images |n|. Here we shall give an alternative 
technique that leads to the result using simple arguments. 

Noting that 

Id d 2 1 / d \ 2 

■^ = 7 > ( 25 ) 



2 dx dx 2 4 \ d\fx 
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the solution of the equation ( j22j ) can be seen as an explicit function of y^i.e. 
R(x,t) = R,(y/x,t), and therefore if R(y/x,t) is a solution, R(— *fx,t) is also a 
solution. Thus, the general solution may be written as a linear combination of 
the two independent solutions 

R(x,t) = aR(s/x,t) + bR(-y/x,t), (26) 

where a and b are determined by the boundary condition at x — as well as by 
the initial condition lim^o R(x,t). 

Reporting (|2^) into (|24|) and using (|l5| ) with ( p6| ) it is easy to show that 
R(x,t) corresponding to the conditions ( |23[ ) reads 



R(x,t) = 



(ex 



eyjn\e(t)\ Jo 
1 



+6 exp 



dz < a exp 



z + 



W)\ V 

Making the change of variable y — ez 2 , (p7|) reads 




R(x,t) 



+ OC 



dyG(x,t;y,0)R(y,0), 



(27) 



(28) 



where G(x, t; y, 0) is the transition probability density given by 



G{x,t;y,0) 



0(ex)9(ey) 
+&exp 



a exp 



(29) 



This is the general solution of Eq. (|22|) satisfying the conditions (|23|). To 
obtain the physical solution we have to distinguish between the two boundary 
conditions at x = : R(Q, t) = for all t and i?(0, i) 7^ for all t. 

For the first case, it implies that the transition probability density from any 
point y 7^ to x — vanishes, i.e. G(x — 0,t;y 7^ 0, 0) = 0. This means that 
the boundary x = acts as a reflecting barrier. This is satisfied by choosing 
a = —6. Furthermore, from Eq. ( |28| ) it is obvious that G(x,t;y,0) must satisfy 
the initial condition lim t _>o G(x, t; y, 0) = 8(ey)S(x — y). Then a = — b = 1. 

Thus, the physical propagator of Eq. (^2j) with reflecting barrier at the 
boundary x = 0, that we denote G lct (x, t; y, 0), reads 



G Tet (x,t;y,0) 



8(ex)9(ey) 

— exp 

6{ex)6{ey) 



exp 



x + y 
■ e(t) 



■ exp 



sinh 



e(t) J 



(30) 
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For the second case, there is non zero probability for a transition from y ^ to 
x = 0. The boundary x — acts then as an absorbing barrier. Again from the 
initial condition G(x,t;y,0) = d(ey)S(x — y) one gets a = 1. From the formal 
solution ( p4] ) it is obvious that if R(x, 0) is a non zero constant in the domain of 
its definition, then R(x, t) does not evolve, i.e. R(x, t) = R(x, 0) for all t > 0. To 
satisfy this condition, it is obvious from ( p8| ) that G(x, t; y, 0) must satisfy the 
"backward" normalization condition 

/+oo 
dyG(x,t;y,0) = l. (31) 
-oo 

Integrating ( |29| ) over the variable y for a — 1 and using ( |3l| ) leads .a = b = 1. 

Thus, the physical propagator of Eq. (|2^) with absorbing barrier at the 
boundary x = 0, that we denote G ahs (x, t; y, 0), reads 



G abs (M;2/,0) 



e{ex)6{ey) 
2y/ne{t)y 

+ exp 



exp 



-^(Vi-^) 2 



1 

e{ex)6{ey) 



■ exp 



x + y 



it) 



cosh 



V e(t) J 



(32) 



Now, we turn to the whole solution of our problem. From the relation (|T^) , 
it is clear that if clq = a(0) ^ (-B(O) 7^ 0), we have 



P(x, 0) = Q(x, 0) = e a °^i?(a;, 0), 



where according to (14), 



Since 



R(x,0) =P(x-a ,Q), 



d 



(33) 
(34) 



P(x,t) = exp (— a(t)) exp — a(i)x— — I exp a(t)— I R(x,t), (35) 



<9.r 



9 



9x 



we get 

P(a:, t) = exp (-£*(*)) i?(x, t), (36) 

with x= (s+ |g)e~ Q W. 

Substituting Eqs ( |3fj| ) and ( |32"| ) for a; = a; into (|2|) and reporting the result 
into Eq. <M) leads 



P£l(x,t) = 6{ex)- 



-a(t) reoc 



dy_ f+*< fsinh ( 2Vgyg \ f 



ey/^{t)Jo Vy* i cosh I e (t) 
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e~°W f 00 dy Jsinh ( 2^/y\ prcf , 



ey/^(f)J-aoVfc l cosh V e(t) y absVy ' 



(37) 



where y = y + o . 

Obviously we may write 



r+oo 

F&(x, t) = / dyKlt{x, t; y, Q)Pg s (y, 0), 



with the propagator given by 
K&(x t t;y,0) - 9 (ex) 9 (ey ) 



e -«(*) 



2v / 7ry e (*) 



, / -s a i - s e Q(t) f (at + y ) \ /sinh ^ 2V^VVo 

>(sx)9(ey ) exp 1 1 ' 1 



\Ary e(t) \ e(i) / \cosh ^ e (t) 



(38) 



It is obvious that lim^o -Klb^O^ *j 2/> 0) = 9(ey)8 (x — y) as it should be. 
However, the "forward" normalization condition on the transition probability 
density is only satisfied for a process with reflecting boundary condition: 

/+oo 
dxK* ei (x,t;y,0) = l, (39) 
-oo 

for alH > and y in the domain of definition. This is obvious since for reflecting 
barrier the process never stop and K Tet (x, t; y, 0) is a proper probability density. 
For a process with absorbing boundary condition, one has 

= f +OC dxK^(x,t;y,0) 



(t) 



= 2m 



(y^)-l + ^exp(--|), (40) 



for all t and y in the domain of definition, with 0K(x) denotes the normal dis- 
tribution. x abs f^fj) i s * ne probability that the process starting at y does not 
reach the barrier before epoch t. In other words it represents the distribution of 
first-passage times at the boundary flC|| . Obviously for this case, K ahs (x, t; y, 0) 
is a defective transition probability density. 
The characteristic functions 

^abs(p) 01 K rcl (x, t; y, 0) and K ahs (x, t; y, 0) with 

respect to x, defined as 

r+oo 

<p5L0>)= dxexp(ipx)Klll(x,t;y,0), (41) 
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may be straightforwardly evaluated. They can be related to Xabs as 
. ( Wn ,.A 

exp 



<Pab 



1(P)=- 



1 1 - ipe(t)e a M y ' 



Xabs 



(1 - ipe{t)e a W)~ 2 " S \ e W ( : _ V e (*) e< * (t) ) / ' 

(42) 

From this expression it is an easy task to obtain the moments of K^ s (x, t; y, 0) 
by the usual formula 

1 B n 

(* n ) (*) = ^^S(P) 



(43) 

p=0 



and check that their corresponding differential equations actually coincide with 
those obtained directly from the differential equation (|l|) with (0) for K^ s (x, t; y, 0). 

Thus we have derived the propagator and its corresponding characteristic 
function for a generalized FP equation with time dependent coefficients satisfy- 
ing the relation ( pfj| ) such that many special cases can be easily deduced from 
it. For instance, if we set B(t) = 0, the propagator obtained in ref.^J is gener- 
alized for arbitrary sign of E(t) and the solution is given for the two boundary 
conditions corresponding to a reflecting barrier and an absorbing barrier. 
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